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Abstract 

The aim of this work is the investigation of the statistical properties of local electric fields in an 
ion-electron two component plasmas for coupled conditions. The stochastic fields at a charged or 
at a neutral point in plasmas involve both slow and fast fluctuation characteristics. The statistical 
study of these local fields based on a direct time average is done for the first time. For warm 
and dense plasma conditions, typically N^. ~ lO^^cm^^, Tg ~ leV , well controlled molecular 
dynamics (MD) simulations of neutral hydrogen, protons and electrons have been carried out. 
Relying on these ah initio MD calculations this work focuses on an analysis of the concepts of 
statistically independent slow and fast local field components, based on the consideration of a time 
averaged electric field. Large differences are found between the results of these MD simulations 
and corresponding standard results based on static screened fields. The effects discussed are of 
importance for physical phenomena connected with stochastic electric field fluctuations, e.g., for 
spectral line broadening in dense plasmas. 

PACS numbers: 52.65.-y, 52.25.Ya 
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The study of the local stochastic electric fields at neutral or charged points in an homo- 
geneous infinite plasma is of interest for several domains. For instance they can be used as 
external perturbations within semiclassical models designed to synthesize line spectra for 
diagnostic purposes 

mm 

. A more general interest comes with the study of non linear 
dynamics of charges undergoing the corresponding forces. This work addresses the simpler 
problem of an hydrogen plasma and the analysis of the local field at neutral points. An 
equivalent study of the local field measured at charged points is straightforward. 

The mass ratio between electrons and protons results in fields involving both fast and 
slow fiuctuations characteristics. In order to analyze the statistical properties of these fields, 
a couple of components, a slow-fiuctuation component (S) and a fast-fiuctuation component 
(F) is introduced, 

E{t) = Ee(t) + Ei{t) = EsloUt) + Efastit), (1) 

where Eg and Ej are the fields due to the electrons and ions, respectively, at an arbitrary 
neutral point. Our objective is to suggest a natural definition of the slow and fast components 
and then determine their statistical properties, subject to the constraints 1) their sum must 
be the total field, 2) they must be statistically independent. Years ago, S and F components 
have been defined in well known articles js, [g], Q, ^ (BMHH) and this will be referred to 
as the standard model. The present work exploits plasma molecular dynamics simulations 
(MD) with a more precise definition that will be referenced as ab initio method. 

The long history of plasma spectroscopy has provided indirect interpretations of the static 
and dynamic properties of the local field through their action on the relaxation of a line 
emitted by the plasma [9|, [lyj. However, the plasma and atomic state averages required do 
not allow detailed information relevant to the splitting of the field into S and F components. 
In contrast, MD simulation provides a means for direct access to such information. The 
results presented here are based on the analysis of the stochastic fiuctuations of sums and 
differences of the ionic and electronic local fields, Ej(t), Ee(t) and the time averaged electron 
field 

1 /-At/a 

^e{t)^, = — E,{t-t')dt' (2) 

l\t J-At/2 



The mean electronic field, Ee(t)^j, calculated on variable periods of time. At, can be con- 
sidered also as a simple measurement device, with variable response time, averaging out on 
the fast fiuctuation of the total field due to the electrons. The reason for introducing this 
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average field is clearly to filter the total field and identify the S and F components Accord- 
ingly, the slow component is defined here as the slowly varying ion field plus the residual 
slowly varying mean electron field 



EsMt) = ^^it)+^e{t)^t■ (3) 



The fast component is the remainder Ei?^At(^) = E(t) — Es,Ati^) — Ee(t) — Ee(t)^j. The 
justification for this decomposition is provided by the simulation results below. 

MD simulation of a partially ionized hydrogen plasma gives the local ion, electron and 
total fields at neutral hydrogen atoms. The density and temperature conditions are chosen 
to explore rather coupled conditions: A^e ~ 10^^ cm~^ and Te ~ 1 eV. Simulations of 
infinite systems are achieved with periodic boundary conditions. Coulomb interactions are 
shielded at a distance A ~ s/2, of the order of the cubic MD cell size s, large enough 
to allow natural screening, e.g., Debye screening, and to not affect appreciably any of the 
properties investigated below. The ion-ion and electron-electron repulsive interactions are 
given by a shielded Coulomb potential Vi2{r) = e^^~~''l^ jr. Attractive electron-ion Coulomb 
interactions are regularized as follows: ViJ^r') = — e^(l — e~'^^^)e~^^'^/r, the short range 
regularization parameter 6 is such that the potential energy of an electron sited on top of a 



proton is equal to the ionization energy of a hydrogen atom HI]. Regularization allowing to 
implement a Coulomb attractive potential in MD simulations has proven useful 12. Il3[|. 

Parameters of interest are the electron-electron or ion-ion average distance, vq = 
(S/AtiN^Y^^, defined in terms of the electron density A"e, the mean electric field modu- 
lus Eq = c/tq, the electron thermal velocity Vq = {kBT^./m^.y^'^, the electron and proton 
coupling constant T = /r^kBTe and the Debye length \o = [ksTe/ A-nNee^Y^"^ . 

Two plasma conditions have been considered, which correspond to a = 0.8 and a = 
0.4, where a = tq/Xd is the parameter introduced in BMHH. Both cases, with a same 
temperature of leV and coupling parameters F ^ 0.2 and F ^ 0.06, correspond to realistic 
experimental plasma conditions. The number of electrons is ~ 1000 with box sizes s > SA^. 
The same time-step t = 1.5 x 10^^ ''s has been used for ions and electrons. This implies 
that, for a = 0.8, 11600 time-steps are required for a proton to move across tq (270 for an 
electron). These parameters give rise to intensive simulations that typically run over a few 
10'' time-steps. Times for ions and electrons to cross the mean distance provide the ionic 
and electronic characteristic times Tj ~ 10~^^s and Tg ^ 10~^^s of fluctuations of local fields, 
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FIG. 1: Field correlation functions for a = 0.8: (Ee(O).Ee(t)), (circles); {^FAt{'^)-'^FAt{t)) for 
At = 0.4Te (crosses); exponential fit (dash line); (Ej(O).Ee(t)), (black circles). 

since the field correlation is lost, in average, once particles have moved across the mean 
distance. A careful preparation of positions and velocities of charges inside the simulation 
cell guarantees a good stationarity of the total energy all along a simulation. Typically, the 
relative variation of the total energy over a simulation is smaller than 1%. The analysis of 
the time averaged electron field introduces new complications beyond those of current two 
component MD simulations and these have been addressed with care. 

In FigH] three field correlation functions have been plotted on the same graph in units 
of Te, for the case a = 0.8: (Ee(O).Ee(t)) (circles), (Ej(O).Ee(t)) (black circles) and 
(EF,At(0).Ej7^At(^)) (crosses) for At = 0.4re. The symbol () denotes a statistical average 
on a relevant set of independent field histories built, according to the present simulation 
method, both at different times and for different neutral points. Clearly, in the framework 
of two component plasmas, the electron field appears inappropriate to represent the F com- 



ponent as the corresponding field correlation 
due to electron-ion coupling mechanisms 10, 



unction manifestly shows a slow de-correlation 



13| , implying a statistical dependency between 



ion and electron fields. In contrast, the correlation of EpAti^) with At = OAt^ is lost over 
a time < Te (estimated only roughly because of the noise). The definitions for the S and 
F components depend on At. It will be understood hereafter that the fast and slow char- 
acteristics of fields can get a more precise meaning by comparison to an additional time r 
connected to some physical process, e.g., the relaxation of a plasma density fluctuation or 
the relaxation of an atomic radiation due to atom emitters imbedded in the plasma. The 
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present S and F component definitions make sense mainly for processes with < t < Ti, 
i.e., for components with fluctuation times roughly connected to those of ions and electrons. 
For instance, cases such that -C r, r <^ re or = oo are not really considered in this 
work. 

The attractive interaction between ions and electrons gives rise to an average response of 
the electrons to a ionic field polarizing locally the plasma. This average response corresponds 
to a natural anti-correlation mechanism i.e. (Ej(O).Ee(t)) < plotted in Fig{T] and also 
clearly shown in [l^. This gives rise to the well-known model based on screening of ion- ion 
interactions by the fast moving electrons. In the standard model the slow component is 
taken to be a screened ion field, with a screening length due to the electrons. A primary 
observation here is that this definition of the slow component does not agree with that of 
Eq.([3]). More precisely, the static screening effect is not equivalent to the time average of 
the electronic motion. A further approximation in applications of the standard model is to 
replace the corresponding fast component by that due to the electrons alone. The results 
below provide significant motivation for reconsideration of those applications, particularly 
in light of the early identification of screening with time averages. 

In order to clarify the following comparisons it is useful to recall a few preliminary notions. 
The Holtsmark function [l^ gives the field distribution function (FDF) inside an infinite 
space filled with a uniform random distribution of point charges Z = ±1 at a density N. 
For instance the Holtsmark curve is the FDF inside an infinite system of noninteracting 
electrons at a density = N while the total ion plus electron FDF, still for noninteracting 
charges, is given by the Holtsmark distribution for the density 2N. It is well known that 
when the interaction between charges are switched on the average field modulus decreases. 
This effect is accentuated by increasing plasma coupling conditions. 

In the two component plasma investigated here, with all the interactions accounted 
for, the electronic and the ionic FDFs are also remarkable. Due to the symmetry of 
charges they are the same. In addition, when At increases, the FDFs of both Ep^Atit) and 
^SAti^) tend towards the common electronic and ionic FDF (lEFDF) as, due to ergodicity 
limAt^oo Ee(t)^^ = (Ee(t)) = 0. The statistical independency of the S and F components is 
effective as demonstrated in Fig{Tl 

Figure [2] is a typical illustration of the central problem of this work. Three FDFs of the S 
component (solid lines), i.e., of the field E5 At(^) = Ej(t) + Ee()f:)^j have been obtained with 
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FIG. 2: S component field distribution functions at neutral emitters, a = 0.8 

1,5-rTT, 1 




|E|/Eo 

FIG. 3: F component field distribution functions for a = 0.8, At increases from the left to the 
right; dash: electron FDF 

various integration times At. The S component shows a small dependence on At around 
Te, suggesting that the smallest time interval contains already most of the contribution 
of electrons to the S component. The strong field wings are located in between the two 
Holtsmark curves (circles) as the slow fields result from interacting charged particles at a 
density 2iVe. The limit of the S component for increasing At is the ion field distribution 
function (dash). As commented above, this ionic FDF resulting from interacting ions at 
a density A^e is shifted towards smaller field modulus with respect to the Holtsmark curve 
(black circles). Finally, for comparison, the field distribution function labelled "Hooper" 
(dotted curve) which represent the standard model, obtained with ionic field at emitters 
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screened at the Debye length, has been plotted. The large difference observed with the 
present model gauges the consequences of the definition that funds the standard model. 

In the same way, the F field distribution function can be extracted from MD simulations. 
In Figj3] when At increases the PDF is shifted from small fields to the electron PDF, i.e., 
the limit obtained when At increases. According to the standard model, the electron PDP 
gives the P component. 

An equivalent behavior is found for the case a = 0.4. Figure H] shows the S field distri- 
bution function obtained with a time-average length At = together with the Holtsmark 
distributions and the ion field distribution. As in the a = 0.8 case the strong field wing of 
the S field distribution function lays in between the two Holtsmark curves. These calcula- 
tions for two distinct coupling conditions show that as expected, the Hooper S component, 
the ionic FDF and the Hotsmark FDF get closer as F decreases. 

Coming back to the scope of this work a few points should be kept in mind. This 
preliminary work exploits intensive relevant MD sampling of stochastic local-fields at neutral 
points in an hydrogen plasma for two distinct density-temperature conditions. Studies for 
other conditions would be useful. However, our objective here, i.e., the study of the statistical 
properties of these sampled field in order to be able to discuss existing models, has been 
reached. Due to the masses involved, local fields undergo an intricate superposition of fast 
and slow fluctuations. This suggests a rational method based on the mean electron field for 
the splitting of the total field into a couple of statistically independent S and F components 
whose statistical properties are investigated separately. The procedure leading to the S/F 
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splitting is straightforward as it basically averages out the fast fluctuation on variable periods 
of time. In contrast, the standard BMHH procedure uses an S component rendered with an 
average ionic potential involving a screening term to account for ions dressed by electrons. 
The S component standard model therefore is a purely ionic field, while that for the ab initio 
method entangles both ions and those of the electrons contributing to slow fields. Clearly, 
the MD ab initio approach is able to show the anti-correlation mechanisms understood in 
the standard model by the screening term implemented in the average ionic potential but 
does not lead to the standard model itself, even as a limit process. The consequence is that 
both approaches result in distinct models for the FDFs. 

Simple comparisons can be performed using the common lEFDF plotted in Figs|2] and E] 
with dash lines. This curve i.e., the limit reached by Es^At{t) and E^? Ai(^) when At increases, 
represents a S and a F component which are not statistically independent as inferred from 
Fig. [H For any At > 0, the S component field distributions are shifted towards stronger 
fields with regard to the lEFDF. In contrast, the F component FDFs are shifted towards 
smaller fields. The same occurs for the standard model whose most probable slow field is 
decreased due to the Debye screening term. In summary, in average F fields are stronger 
than S fields for the standard model. The opposite is found with ab initio methods implying 
that the balance between the S and F fields results drastically different for the present and 
the standard models. The unexpectedly large magnitude of the observed differences of the 
FDFs should be interpreted as a warning in favor of a careful use of the S and F components 
concepts, for instance, in the context of Stark broadening of lines. 
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